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We propose a general scheme for the derivation of the signals resonant inelastic (and elastic) x-ray 
scattering (RIXS) gives access to. In particular, we find that RIXS should allow to directly detect 
many hidden orders, such as spin nematic, bond nematic, vector and scalar spin chiralities. To do 
so, we choose to take the point of view of effective operators, leaving microscopic details unspecified, 
but still keeping experimentally-controllable parameters explicit, like the incoming and outgoing po¬ 
larizations of the x-rays. We ask not what microscopic processes can lead to a specific outcome, but, 
rather, what couplings are in principle possible. This approach allows to systematically enumerate 
all possible origins of contributions to a general RIXS signal. Although we mainly focus on magnetic 
insulators, for which we give a number of examples, our analysis carries over to systems with charge 
and other degrees of freedom, which we briefly address. We hope this work will help guide theorists 
and experimentalists alike in the design and interpretation of RIXS experiments. 


I. INTRODUCTION 


Many systems have ground states with well-defined 
order parameters which couple directly to conventional 
probes such as neutrons or light. The accessible data 
usually comes in the form of “structure factors,” i.e. cor¬ 
relation functions of two “elementary” observables. Clas¬ 
sic examples are magnetically ordered states, e.g. ferro- 
magnets and antiferromagnets whose magnetic structure 
and fluctuations can be resolved by methods like neutron 
scattering, muon spin resonance (/rSR), nuclear magnetic 
resonance (NMR) etc.. However, many of the “exotic” 
phases proposed by theorists do not fall into this cate¬ 
gory. Some states exist, for example, which possess a 
well-defined local order parameter, but still evade robust 
characterization using “conventional” probes. The order 
is then commonly referred to as “hidden.” Typically, the 
order parameters of such systems have quantum numbers 
which are multiples of those which elementary particles 
give access to when coupled linearly to the system. For 
example neutrons can excite S = 1 magnons, but not 
S = 2 excitations (owing to the dipolar coupling between 
the neutron and electron’s spins). Perhaps the simplest 
and best-known example of a hidden order is that of spin 
quadrupolar (also called nematic) order [T]. In that case, 
the expectation values of the spin projections, (S'f) (note 
the spins transform as “dipoles”) are zero, but those of 
“quadrupolar” operators, like (S^S^), are not. Many 
other types of hidden orders have been proposed in the 
literature. Among those are spin “bond nematic” 00 , 
where the order parameter contains spins on neighboring 
sites, and spin vector and scalar chiralities, which involve 
antisymmetric products of spins. Hidden orders also arise 
in conducting systems, with the famous example of ne¬ 
matic (in that case, “nematic” refers to rotation -discrete 
or continuous- symmetry breaking in real space) order in 
the pnictide superconductors. 

Here we show that Resonant Elastic and Inelastic X- 
Ray Scattering (REXS and RIXS) can in principle mea¬ 
sure spin nematic, vector and scalar chirality, and many 
more correlation functions (static and dynamical for 


REXS and RIXS, respectively). In general, we propose 
an enveloping scheme which allows to systematically enu¬ 
merate which correlation functions will contribute to the 
RIXS signal in any given polarization geometry. REXS 
signals are obtained from RIXS in the w —> 0 limit. In 
particular, in the case of static order, REXS signal should 
display corresponding “Bragg” peaks. 

“Resonant scattering” refers to techniques where the 
energy of an incoming probe is tuned to a “resonance” 
(a.k.a. “edge”). In that case, not only is the absorp¬ 
tion (virtual or real) cross-section dramatically enhanced, 
but the latter may also involve nontrivial operators, al¬ 
lowing to probe correlation functions of complex order 
parameters, i.e. typically those of hidden orders, which 
are otherwise hardly accessible. This is clear upon think¬ 
ing in terms of perturbation theory in the probe-system 
coupling amplitude, and we soon specialize to an x-ray 
probe. The scattering amplitude up to second order is 
given by 00 


r fi = <f|tf'|i> + E 


(f\H'\n){n\H'\ 
E± - E n 


(1) 


where |i,f) denote the initial and final states of the 
{system + electromagnetic (EM) field}, H' is the cou¬ 
pling Hamiltonian between matter and the EM field, 
{|n)} forms a complete set of states (the “important” 
ones will be discussed later) of the system, and E a is the 
energy of the state |a). When there exist states | n) which 
are close in energy to E ±, the system is said to be at reso¬ 
nance with the probe and the second order amplitude in 
Eq. (|T|) largely dominates the first. Moreover, within per¬ 
turbation theory, the former contains, among others, the 
following chain of (virtual) processes: the absorption of 
a photon, the evolution of the resulting system, followed 
by the emission of a photon. The RIXS signal is the 
cross-section relative to the amplitude of such a process, 
when the incoming x-ray light is tuned to a resonance 
which involves the excitation of a core electron to a va¬ 
lence level, i.e. when | n) is a state of the pure system (no 
photons) and contains a “core hole”. Typical orders of 
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magnitude for such x-ray frequencies range between 0.01 
and 10 keV [SHZ], he. correspond to photon wavevectors 
of order 1-10 -3 A -1 . 

Detailed microscopic analyses of RIXS processes in a 
number of systems have been described at length in the 
literature [5], some even predicting the observation of 
correlation functions of complex order parameters mm- 
Here we do not belabor on them, but rather base the 
analysis entirely on the observation that the initial (be¬ 
fore the photon is absorbed) and final states (after the 
photon is emitted) of the system both belong to its low- 
energy manifold. Essentially, in that approach, the only 
important feature of the microscopies is the reduction of 
(at least spatial) symmetries to those of the core-hole site 
point group. Such a symmetry-based strategy has a few 
major advantages. An accurate description of all possi¬ 
ble microscopic processes is a very complex many-body 
problem, which is moreover subjected to many uncertain¬ 
ties concerning the atomic structure in a material. As 
a consequence, such approaches are inherently material- 
specific. It is moreover very difficult to exhaust all pos¬ 
sible processes through microscopic reasoning. The sym¬ 
metry procedure bypasses these issues. This type of fully 
effective approach was recently insightfully pioneered in 
Ref. GDI in the context of magnetic insulators, where the 
author gave the form of on-site effective RIXS operators 
for up to two on-site spin operators. Here we construc¬ 
tively rederive and generalize Ref. [T0}s main result to 
all possible symmetry-allowed couplings, including those 
which involve multiple-site operators and degrees of free¬ 
dom other than just spins. Moreover, the broader con¬ 
text of the derivation presented here helps make more 
transparent the correlations possibly probed in RIXS, on 
which we focus. 

The remainder of this paper is organized as follows. We 
first review the form of the light-matter interaction, the 
important symmetries to be considered, and derive the 
form of the effective operators whose correlations RIXS 
measures in insulating magnetic systems, which are sum¬ 
marized in Table[l] We then turn to the study of three im¬ 
portant examples of hidden orders as may be realized in 
real materials: spin nematic order, bond nematic order, 
vector and scalar chiralities, and calculate the expected 
RIXS signals in these three concrete cases. At the end of 
the paper we briefly address systems with charge degrees 
of freedom. 


II. EFFECTIVE OPERATORS 


The leading order Hamiltonian H' which couples light 
to matter and is involved in the second-order ampli¬ 
tude of the interaction cross-section is given by, in the 
Coulomb gauge V • A — 0 •"> [TT] : 


H' = y: — -$ r • A r + ipl xj) r • (v X A r 

' m 2m \ 


( 2 ) 


with the vector potential 


E 


fl V "A 

2VeoCJk ^^ 


*^t ikr 


£ a } 


k,e c 


+ h.c. 


( 3 ) 


H' acts in the product space of the electrons TL e - and 
photons Hphot , W = H e - x Hphot, and xf’ are the 
electron creation and annihilation second-quantized op¬ 
erator fields, h is Planck’s constant over 27 t, e and m 
are the electron charge and mass, respectively, a* and a 
are the photon creation and annihilation operators, e de¬ 
notes the photon polarization, V is the volume in which 
the EM field is enclosed, eo is the dielectric polarization 
of vacuum and Wk = k = c|k| where c is the speed 
of light. Here, for concreteness, we make two approxi¬ 
mations, often used in the literature [5j: we consider (i) 
that |k• <5r| <C 1 at the relevant x-ray wavelengths, where 
r = R + dr where R denotes the position of a lattice site, 
and so, at zeroth order, e tk Sr ~ 1 [I2j, and (ii) that in 
Eq. § the magnetic term (oc er) is subdominant com¬ 
pared to the “electric” one (oc p). We return to these 
approximations in Appendix [E] Therefore, the second- 
order RIXS amplitude for processes with a core hole at 
site R reduces to 


7S'= X < f 
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where O 



pGp, with G = 


Ei+hojq—E n 5 


where | n) are restricted to “intermediate” states with a 
core hole at site R (i.e. close to resonance). The second 
expression Eq. <© is obtained by requiring |i) = |z) ( 8 ) |ker) 
and |f) = |/)(g> |k'£ 7 ). Importantly, O acts purely in elec¬ 
tronic space, and moreover within the low-energy man¬ 
ifold, provided the system immediately “returns” to a 
low-energy state as the outgoing photon is emitted, as is 
usually assumed. We therefore ask: what effective op¬ 
erator acts purely in this manifold which reproduces the 
matrix elements 7^ f ? If we know the low-energy mani¬ 
fold and a basis which spans it, and if the basis elements 
are physically meaningful, we shall immediately obtain 
which correlation functions RIXS produces. We insist 
once again that, within this approach, all “intermediate 
processes,” no matter how complicated, are in a sense 
included, and need not be discussed. 

As usual, most general arguments stem from symmetry 
considerations, which we now address. The core hole is 
immobile, which imposes a strong symmetry constraint 
on O r : it should be invariant in real space under point 
(site R) group symmetries. Another constraint comes 
from the “locality” of the effect of the core-hole in the 
“intermediate propagation time” r = 1/T ~ 10“ 15 s [5], 
which implies that only operators which act in close prox¬ 
imity to the site of the core hole should be involved. 
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While this statement may appear somewhat loose, a 
quick order-of-magnitude analysis shows that, even in 
a metal , electrons will not travel over more than very 
few lattice spacings over the time r [T3L Finally, since 
transition amplitudes are scalars, by keeping the polar¬ 
ization dependence explicit, we impose constraints on the 
combination of operators which multiply the polarization 
components. This is what we address now and is sum¬ 
marized in Table HI 

For concreteness and ease of presentation of the deriva¬ 
tion we now focus on magnetic insulators, though we note 
that the same ideas carry over to systems with charge 
(and other) degrees of freedom, to which we return at 
the end of the paper, in Sec. |Vj Indeed, because of the 
“locality” of the effective scattering operators, insulating 
systems are more readily tackled. Local (in the sense 
of acting only on degrees of freedom living in a small 
neighborhood in real space) operators in insulating sys¬ 
tems yield a very natural description of the system, and 
the low-energy manifold, being finite (generally a well- 
defined J multiplet, possibly split by crystal fields) and 
sharply defined (usually a gap separates multiplets), can 
be spanned by effective “spin” operators (finite vector 
spaces of identical dimensions are isomorphic). Therefore 
only a spin operator basis compatible with the combina¬ 
tions of polarizations remains to be found. 

In the absence of both spin-orbit coupling at low ener¬ 
gies (core levels always experience very strong spin-orbit 
coupling i) and of a magnetic field, the system should 
be rotationally symmetric in spin space. Moreover, in 
principle, in the Hamiltonian, under spatial symmetries, 
the spins are left invariant. However, here, in the RIXS 
structure factor, the situation is more subtle. Spin ex¬ 
citations (and hence spin operators) may only arise in 
the structure factor thanks to spin-orbit coupling at the 
core. Therefore, in principle the structure factor itself 
should display signs spin-orbit coupling urn nn, with 
the effective spin operators transforming under lattice 
symmetries. Even upon neglecting transition operators 
which break rotational symmetry if spin orbit coupling is 
weak at the valence level, the effective spins still trans¬ 
form under real space symmetry operations. .175] Then, 
the polarizations and (effective) spins (the latter make 
up the operators 0 Miy , as mentioned above) transform 
as usual vectors and pseudo-vectors, respectively, under 
spatial transformations, and according to e —> —e* and 
S — S under time reversal (see Appendix [B]). In other 
words, under the full spherical symmetry group, using 
the notations from Ref. □1 e and S transform under 
D and D +, respectively (under SO( 3) operations, both 
the polarization and spin vectors transform under the 
L = 1 representation). Since Dp x Df = + D+ + D (j~ 

(lxl = 0 + 1 + 2 for SO(3)), any combination of spin op¬ 
erators which transform under the same representations 
can in principle be involved in the RIXS signal. Depend¬ 
ing on the number of neighboring operators one chooses 
to include (and on the value of S(S + 1)), possibilities 
differ. The situation for up to three spin operators (on 


the same or nearby sites, from “locality”) is summarized 
in Table [I] (see in particular the caption), and details of 
the derivation are given in Appendix [C] 

On-site terms. — Upon considering on-site terms only 
(i = j = k), where one need not take into account any 
further lattice symmetries, and up to two spin operators, 
we recover the expression from Ref. urn 

Ti = a 0 (s'* ■ e) + ai (s'* x e) ■ S, + a 2 [ e '*, e] [S*, S t ], (6) 

where T) = £'*0^ v £ v , and where [Sj,Sj] is the trace¬ 
less symmetric second rank tensor constructed from S, 
and S j, i.e. given by: = \ (SfS'J + S”S%) - 

!(Si-Sj)<5 M! ,, and analogously for [e 7 *, e]. The symmetric 
product [£'*,£][Si, Si] = has 

all indices contracted. The a n are material-specific co¬ 
efficients pH] . The generalization to discrete symmetries 
is formally straightforward (though usually gruesome in 
practice) and discussed in detail in Appendix |D| 

Off-site terms.— The above considerations take care 
of the symmetry aspects relative to spin space. To fulfill 
the constraints associated with the lattice, which enters 
through S r —>• [det R) R •s*. r where R is a spatial op¬ 
eration (see Appendix the expressions must be ap¬ 
propriately symmetrized. For example, take a Id chain 
of S = 1/2, and consider a maximum of two spin terms. 
Then, if lattice sites are centers of inversion, the transi¬ 
tion operator will be (still assuming spherical symmetry 
in spin space): 

Ti = a 0 (e'* • e)S, • (Sj_i + Sj+i) 

+( £/ * x £ )' ( a i,iSi + cci^Si x (Sj_i + S i+1 )) 
-\- 0 i 2 \_e'*, e][Sj, Sj_i + Sj+i], (7) 

where the a n and a n ^ m are material-specific coefficients 
which multiply terms which belong to the same irre¬ 
ducible representation (n) (or copy (to) thereof if an ir¬ 
reducible representation appears multiple times). 

From Table |TJ one may directly read out the quantities 
whose correlation functions will contribute to the RIXS 
signal, as well as which polarization geometry will reveal 
them while switching off (most of) the other contribu¬ 
tions (e.g. e'* || e will “switch off” the s'* x e “channel”). 
Indeed the differential cross-section is given by |4] 
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where 50 and 5E denote elementary solid angle (related 
to the momentum transfer k —k') and energy, respec¬ 
tively, and where A E is the measured energy transfer. 

Before moving on to the discussion of specific exam¬ 
ples, we make a couple of important remarks, (i) It is 
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representation 

polarizations 

one spin 

two spins 

three spins 

0 

e f * • e 


Si • S j 

(Si x Sj) ■ S k 

1 

e'* x e 

s* 

Si x Sj 

(Si -Sj)Sfc, (Si x Sj) x Sk, [Si, Sj] ■ St, 

2 

[s'*. ®1 


I&.Sj-l 

[Si x Sj, Sk], [Si, Sj] x Si, 


TABLE I. Generic form of the operators in magnetic systems which couple to combinations of the polarizations in the absence of 
spin-orbit coupling. The double brackets represent the traceless symmetric products, [u, v] M!/ = | + u v u M ) — | (u ■ v)S flly , 

and the dot and vector products between a matrix (obtained through [.,.]) and a vector are defined such that: ([u, vj • 
w )m = XZ„I U , vlpvWv, ([u, v] x w) w = Yjv,\ e^A P [[u, (see Appendix 0). Moreover, the product [u, v][w, t] also 

denoted [u, v] ■ flw, t] is defined to be the fully symmetric product with all indices contracted: „[u, v] M „[w, t] M „. Each 

row corresponds to a given irreducible representation of a combination of the incoming and outgoing polarization, given in 
the second column. Each entry on the right of the double bar gives the combinations of spins which transform as does the 
combination of polarizations on the same line. The columns simply indicate the number of spin operators involved in the 
effective operator. In principle, RIXS may measure correlations functions of the operators given in this table. This table is also 
“valid” for matrices which connect local “band” indices with the same symmetries in conducting systems. See Sec. E3 


important to note that, for effective spin-1/2 systems, 
only off-site terms can contribute to, for example, the 
[e , *,e] channel. Indeed, there exist only four (count¬ 
ing the identity) linearly independent S' = 1/2 opera¬ 
tors. Therefore, while off-site contributions are expected 
to be weaker (they may only arise from so called “indi¬ 
rect” processes 0 ), in an effective S = 1/2 system, a 
“multi-site” signal in the [e'*,e] channel will not “com¬ 
pete” with signal from possibly-larger onsite couplings, 
offering hope to unambiguously detect such correlations. 
(ii) We caution that, of course, this symmetry-based ap¬ 
proach does not any give information on the absolute or 
relative strengths of the signals in the different channels. 
Moreover, “selection rules” relative to the chosen “edge” 
need to be additionally taken into account, (in) An ad¬ 
ditional word of caution is in order: as far as we under¬ 
stand, the measurement of the outgoing polarization is 
not currently possible on instruments being used at this 
point, although the new state-of-the art facility currently 
under construction (which will also provide much higher 
resolution in energy, currently at around 100 meV) will 
be able to. 


III. SPIN NEMATIC IN THE 
BILINEAR-BIQUADRATIC S = 1 MODEL ON 
THE TRIANGULAR LATTICE 

The 5=1 bilinear biquadratic model with Hamilto¬ 
nian 

J? = ^(J 1 S i -S, + J 2 (S l -S,) 2 ), (9) 

(id) 

on the triangular lattice has been quite intensively stud¬ 
ied, especially so in recent years after it was suggested 
that it could be relevant to the insulating material 
NiGa 2 S 4 , where Ni 2+ is magnetic, with S = 1 [UK 
[23] . This material is made of stacked triangular planes 
of Ni 2+ ions, and displays no long-range spin ordering 
but low-temperature specific heat which grows with tem¬ 
perature as T 2 [18] , The latter facts motivated the 


minimal description of NiGa 2 S 4 by the model Eq. ©>• 
which, for Ji > 0, features two quadrupolar phases, 
one “ferroquadrupolar” and one “antiferroquadrupolar.” 
These phases are characterized by a vanishing expec¬ 
tation value for the spins, (5/) = 0, but an on-site 
“quadrupolar” (a.k.a. “spin nematic”) order parameter: 
({5/, 5/} — 2^ 0 (a diagonal part is subtracted 
to obtain a traceless operator). Since here we look 
not to accurately make predictions for the actual mate¬ 
rial NiGa 2 S 4 , but rather to demonstrate that RIXS will 
provide unambiguous signatures of quadrupolar order, 
we now restrict our attention to the minimal bilinear- 
biquadratic model Eq. © , despite the fact that the latter 
will clearly not account for all the experimental features 
(not discussed here) of NiGa 2 S 4 [22] . 

The wavefunctions of nematic states are simple single¬ 
site product wavefunctions. For spin-1 systems, prod¬ 
uct wavefunctions can generally be expressed as | ip) = 
IL d ; ' lu), where we have defined |rj) = (| x { ), \ y*), \z,)), 
where d; £ C 3 and |d, | = 1. The states | /r,;) are time- 
reversal invariant states defined such that 5f |/z») = 0, 
i.e., in terms of the usual eigenstates of the Sf operator, 

l*> = ^(|l> - IT)), I y) = ^(|i> + l T » and I*) = -*l°> 

[Mj . In the case of a “pure” quadrupolar phase, for this 
basis choice (with time-reversal invariant states), d^ £ R 3 
[24] , which one can check indeed leads to (Si) = 0. The 
vector di at each site is called the “director,” and corre¬ 
sponds to the direction along which the spins do not fluc¬ 
tuate. In nematic states, the direction along which the 
director points may vary at each site, like in the “antifer¬ 
roquadrupolar” phase of the above model, where the di¬ 
rectors form a three-sublattice 120° configuration. In the 
ferroquadrupolar phase, the directors on each site point 
in the same direction, which can be arbitrarily (since the 
Hamiltonian is isotropic in spin space) taken to be the 
z direction. In that case, the unit cell is not enlarged. 
In ordered (or field-polarized) ferromagnets and antifer- 
romagnets, the low-energy elementary excitations of the 
system are spin flips/waves, i.e. S z = ±1 local excita¬ 
tions. In nematic states, where it is the directors which 
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are ordered, spin waves translate to “flavor waves” where 
there are now two pairs of conjugate “transverse” bosons. 
Flavor wave spectra and dipolar and quadrupolar corre¬ 
lations for the model Eq. © on the triangular lattice 
have been calculated in several references n [ini ns on. 
Our derivation is provided in Appendix |F 1[ and here we 
give the full RIXS structure factor for the model, assum¬ 
ing on-site spin operators only (expected to provide the 
largest contributions to the signal), and spherical sym¬ 
metries (a derivation is provided in Appendix |P|), and 
provide a few plots in Figure [T] for various polarization 
geometries and assumptions on relative absorption coef¬ 
ficients (about which symmetry analysis gives no further 
information). 


Z RIXS oc 
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is the second rank fully 
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where A q = 

7/(^2 — Ji ), Wq 

2 (cosgi + cos(|[gi + %/3g 2 ]) - 
and = aie ll x p £ x £ , * p (e 

(O') 

antisymmetric tensor), = 02 (—2/3 

£ ,p *£ v -\-£ p ‘*£' 1 ') (note that aq and a 2 depend, in partic¬ 
ular, on the details of the atomic and crystal structures 
see Appendix |F 1[ Quadrupolar correlations are 
therefore directly seen. Clearly, one recovers the proper 
scaling of the amplitudes for the Goldstone mode (the 
system spontaneously breaks spin-rotation symmetry 
in the ferroquadrupolar phase) at q = 0 at low energy, 
w q - |q| and J RIXS - f erro ^ i/ Wq . Figure [I] illustrates the 
associated smoking gun evidence for quadrupolar order 
provided by RIXS. 


IV. BOND NEMATIC AND VECTOR 
CHIRALITY IN NEAREST AND 
NEXT-NEAREST NEIGHBOR 5=1/2 
HEISENBERG CHAINS IN A FIELD 


The 5 = 1/2 ferromagnetic nearest-neighbor and an¬ 
tiferromagnetic next-nearest-neighbor Heisenberg model 
on a chain 

H = J2 • S ‘+1 + J 2 S * • s *+2 - hS?) (11) 

i 

with J 12 > 0 is thought to be a minimal model for 
LiCuVCH, a distorted “inverse spinel” (with chemical for¬ 
mula ABB’ 04 ) material such that the system can be seen 
in a first approximation as a set of parallel edge-shared 
Cu0 2 chains separated by Li and V atoms . Cu 2+ 

are magnetic ions with spin 1/2. As will be important 
later, we note that the point group symmetry at each 
Cu site contains inversion symmetry. This material dis¬ 
plays a complex phase diagram, which is now believed to 
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FIG. 1. Color plots of the a) spin-spin correlation function 
(as probed by e.g. inelastic neutron scattering or RIXS with 
e 1 * _L e and 02 small enough), and b) signal probed by 
RIXS for s = (i, l,0)/\/2 and s' = (1, z, 0)/x/2, both for the 
model of Eq. © with J 2 /J 1 = — tan(7-7r/16) (ferroquadrupo¬ 
lar phase ESI) on the triangular lattice, c) and d) Equal time 
(integrated over frequency) versions of the signals shown on 
plots a) and b), respectively. Note that the intensities are 
independently normalized. On figure b), the intensity is seen 
to diverge at the Goldstone mode, in sharp contrast with the 
vanishing of the spin-spin correlation function at the same 
point in figure a). In plots a) and b), Cj = oj/^/J\ + Jf. 


show, from low to high field: incommensurate helical or¬ 
der, spin density wave order along the chains, and, possi¬ 
bly, right below the saturation field, a spin nematic state. 
Why the J\ — J 2 Heisenberg model of Eq. seems like a 
reasonable starting point to describe this material may be 
articulated as follows: (i) there is experimental evidence 
for chain structure physics (see above), (ii) Cu usually 
displays weak spin orbit coupling, suppressing any strong 
anisotropy in spin space, and (iii) further-neighbor inter¬ 
actions in such compounds are usually sizable, owing to 
the configurations of the exchange paths. In fact, J\ and 
J 2 were estimated to be 19 K and 44 K, respectively, us¬ 
ing neutron diffraction and susceptibility data on single 
crystals 1271 I30j . Moreover, in some parameter regime, 
a number of the phases numerical simulations obtain for 
the model are reminiscent of those experimentally ob¬ 
served in LiCuVCH, as we now discuss. 

For J 2 /Ji > 1/4, in a non-zero but weak enough field, 
the minimal model has been shown to exhibit a nonzero 
vector spin chirality z • Xi,i+i = z • (Sj x S^+i) and z • 
Xi,i +2 = z • (Sj x S, +2 ) (a non-zero z-component of the 
chirality does not break any continuous symmetry of the 
model in a field applied along the z-axis and is therefore 
allowed), reminiscent of the helical order in LiCuVCH. 
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More precisely, DMRG and exact diagonalization have 
probed signs of long-range chirality correlations |3H - l33l| , 
and the bosonization of the field theory—which unveils 
a Luttinger liquid phase— predicts (x*,i+i • z) y 0 and 
(Xi,i+ 2 • z) y 0 J32H311- The higher-field phase of the 
model numerically shows evidence of (bond) quadrupolar 

correlations HI 131133]. 

Again, here we claim not to provide a detailed descrip¬ 
tion of the material, but we propose that RIXS might 
be able to probe vector chirality as well as bond-nematic 
order in this system. 

In order to compute the RIXS signal, we proceed like 
in Ref. [323 closely follow their derivation, and start from 
the limit J\ <C Ji of two decoupled chains (each with lat¬ 
tice spacing 2ao). Each one may then be independently 
bosonized. We use the conventional notations for the 
boson fields, and <(> 1 , 2 , where the indices are chain 
labels, and \<j) v {x', t'), d x 6fj,(x, r)] = iS^ v 6(x — x')S(t — 
t'), for 5 , v = 1,2. The spin operators are given by 

HEDE2U311 

S+(x) = e l ^ e ^ x) (12) 

x ^(—iyb + b' sin(27r Mj + 4>fj,(x))j 

S z (x) = M + ~^=d x (j)^{x) (13) 

V 71- 

— (—l)-'asin(27rAfj + -\/47r</> M (:r)), 


where x is the coordinate of a site, while j £ Z la¬ 
bels a “unit cell” of two sites { 1 , 2 } (sites can be la¬ 
belled by l = 2 j + n), M is the total magnetization 
(due to the field), and a, 6 , 6 ' are non-universal con¬ 
stants. Note that here the subscript 5 in 5“ is un¬ 
related to the subscript i in Eq. (11). As mentioned 
above, when Ji = 0 , the two chains are decoupled and 
each one obtains two free-boson theories, with the action 
S cs = I dx ! <M| ( K (d x 0u) 2 + ^{9 x (j> IJ .) 2 )+id x e i fAt/v], 
where K and v are the Luttinger liquid parameter 
and spin velocity of the antiferromagnetic Heisenberg 
(J 2 ) spin chain in a field. Ji / 0 introduces cou¬ 
plings between the chains. Then it is useful to define 
7 ± = (71 ± 72 )/v / 2 for 7 = 0,<j). The coupling ac¬ 
tions are 6 >i = <71 f da; sin (v / 87r</>- + nM) and S 2 = 
g 2 J dx{d x 9j r ) sin a/27t 6L, where 0 <M<l/ 2 is the mag¬ 
netization per site, with parameters g\ = — J 1 a 2 sin 7 rAf 
and g 2 = —Jiy/nb 2 /y/2, which will lead to “bond ne¬ 
matic” and “vector chiral” phases. This model displays 
scale invariance, and renormalization group (RG) ideas 
apply. Then, within this approach, if < 72/51 Hows to zero 
(resp. infinity) under the RG flow where high-frequency 
modes are integrated out, the system goes into the ne¬ 
matic, where cf>- gets pinned to a value which minimizes 
the integrand of <S 1; (resp. vector chiral, where it is the 
integrand of S 2 which acquires a finite expectation value) 
phase [32]. Details are given in Appendix F 2 


Because each site on the chain has only two neighbors, 
we expect that the contributions to the RIXS signal from 
three-spin interactions should be extremely weak. So, 


from Table [TJ assuming a weak enough effect of spin-orbit 
coupling at the low-energy level, the RIXS transition op¬ 
erator is given by Eq. 0 in zero held, and by 

Ti = ao.jVx* • e ± )S± ■ (S^i + S/ +1 ) (14) 

+aoAe , z *e z )SnS z _ 1 + SUi) 

+(e'* x e) z Si + a lt 2,zSi x (S,_! + S. l+1 )) z 
+a 2 „ L [e , *,e] J -[S I ,S l _ 1 + S J+1 ]\ 


for h y 0, i.e. when the full SU (2) symmetry in spin space 
is broken down to 17(1) (see Appendix [D]) . In Eq. (141, 
we used the definitions u = u_l + u z z and [u, v]/, = 
\{u^v v + /u' 1 ) — |( u ± • vi)^ with n,u = x,y only. 
The a n ^ and o; n)TO;jLt are coefficients. Finally, we fold 
the following low-energy (long distance and time) leading 
contributions (see Appendix |F 2\ to the RIXS structure 
factor: 


7-nematic 

L co,q 


OC 


E 

e=±l 


0(u ; 2 — v\(q — e7r) 2 ) 


^ Ju 2 -vl(q-eir) 


-2-1/K+ ’ 


(15) 


for, e.g., e x e'* = 0 and e 1 z in the nematic phase, 
and, around q = ±2nM ± n: 


Z^q al oc ^2 B(w 2 - v\(q - 2neM - e'-rr) 2 ) (16) 

e,e'=±l 

/- 1 '' -2 

x/w 2 — v+(q — 2neM — e'n) 2 

in cross polarizations, with (e x e'*) || z in the vector 
chiral phase. In the expressions above, K + = K(1 + 
Ji^) and v + = v{l - Ji^) [note that K(M = 0) = 1/2 
and K{M = 1/2) = 1], Figure [ 2 ] displays some examples. 



q-n q-TT(2M-‘\) 


FIG. 2. Color plots of the dominant contributions to the 
a) [e , *,eJ _L channel (fourth line of Eq. ( |14| )) in the nematic 
phase, around q = n, as given in Eq. (15 1 b) (s'* x e) z channel 
(third line of Eq. ( |14} ) in th e ve ctor-chiral phase, around q = 
7r(2M — 1), as given in Eq. (161. 


V. OTHER DEGREES OF FREEDOM: 
ELECTRONS, PHONONS AND ORBITALS 

The derivation of effective RIXS operators presented 
above in the context of magnetic models readily extends 
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to systems where other degrees of freedom are important. 
Indeed, the symmetry arguments we employed are gen¬ 
eral enough that they carry over to any type of problem. 

Modifications arise at the level of the identification and 
choice of basis for the space of operators which act on the 
local low-energy manifold. In magnetic insulators, as dis¬ 
cussed earlier, the natural degrees of freedom are on-site, 
and a Hamiltonian is always associated with the spec¬ 
ification of what the local degrees of freedom, namely 
effective “spins,” are. More microscopically, one can see 
an effective spin degree of freedom “emerge” from the 
multiplet structure of a single-ion Hamiltonian at each 
site. Now, similarly, if orbital degrees of freedom are to 
be treated explicitly in an insulating system in RIXS, 
one may simply introduce a set of (effective) operators 
L, L^L 1 ' etc., which transform as pseudo-vectors under 
real space operations, and obtain a table similar to Ta¬ 
ble |TJ where now each row should be associated with an 
irreducible representation of the appropriate point group. 

Now, systems with charge degrees of freedom, or 
phonons, are usually approached from a more field- 
theoretic perspective, where one has lost sight of a mi¬ 
croscopic model, and operators are labeled by some mo¬ 
mentum index (among others). That being said, given 
a material, one may always, much like for the insulat¬ 
ing magnet case, think about how many electrons, and 
which single-ion orbital (or spin-orbital), a given ion will 
“contribute/provide” to the valence band of the whole 
solid. Provided one can determine this, it is reason¬ 
able to think of these spin-orbital states and number of 
electrons as the building blocks for the local low-energy 
manifold relevant to RIXS, and the basis of operators 
can be made of those which reshuffle the electrons in 
the (single-electron) spin-orbital states (even if the elec¬ 
trons interact, such a non-eigenstate basis can be cho¬ 
sen nevertheless). As an example, consider an atom 
A contributes n on-site states to the valence band(s) 
of the system, with creation and annihilation operators 
tpl a ,ip ra- One may build on-site operators tf>l a M a pip r p, 
‘•Plai’lpMotP'ySipr'yi’rS etc., where a,/3,7,<5 = 1 ,..,n (may 
be orbital labels, for example), as well as some involv¬ 
ing neighbors, il’l a M a pip r 'p, V’Jaf'Vetc.. 
Despite the more delocalized nature of the electrons in 
an itinerant system, a quick order of magnitude estimate 
shows that, even in a typical metal, only close-neighbor 
operators are involved in the RIXS transition operators 
(see Sec. [0] and foonote therein). An additional con¬ 
straint in RIXS is charge conservation, since no electrons 
are kicked out of the sample. Then, much like in the 
case of magnetic insulators, we may split the tensors M 
into irreducible representations and obtain the coupling 
terms to the corresponding combinations of polarizations. 
In a single band model, for example, the only on-site op¬ 
erators are the density ■i/’JV’r and spin iplcr'i/t r j55] (and 
powers thereof, though the latter should be expected to 
contribute sub-dominantly). 

Like in any endeavor to compare experiment with the¬ 
ory, in any other technique, the most-delicate step in the 


calculation of a structure factor in a given ground state 
will be to understand how the ip TCt operators from the ba¬ 
sis act on this ground state and are related to quasiparti¬ 
cle operators, if any. This is particularly true in the case 
of metals (but also of course in that of, e.g., quantum spin 
liquids), where, even in the case of a Fermi liquid, where 
the notion of quasiparticles is meaningful, the quasipar¬ 
ticle operators are, in the crudest approximation, re¬ 
lated to the electron operators through the square-root 
of the quasiparticle weight 0 < Z < 1: ~ \/~Z^. 

Therefore, a factor of at least Z 2 will be involved in the 
contribution of a quasiparticle-related excitation to the 
RIXS cross-section. Because Z can be very small, like in 
a highly correlated metal, it is important to keep track 
of those factors to estimate the (esp. relative) amplitude 
of a signal of a given origin. For example, upon taking 
the a minima point of view of a single-band Fermi liquid 
[36l for the overdoped cuprates, one should keep in mind 
that factors of Z are likely to greatly suppress the quasi¬ 
particle contribution to the RIXS signal. This should be 
crucial in deciphering the origin of the features seen in 
RIXS spectra of those materials (3171 - 711)1 . 

The case of phonons is quite similar. At the symme¬ 
try level, phonons bear no spin degree of freedom, but are 
associated to lattice degrees of freedom and their symme¬ 
tries. There may be several phonon/displacement modes 
at each site, so that one can introduce several phonon 
creation operators cj a . The symmetries to be considered 
should be purely spatial, and related to point group sym¬ 
metries at site r. Phonons and orbital degrees of freedom 
are likely to be important in the context of the nematic 
order seem in the pnictide superconductors, whose mi¬ 
croscopic origin is not yet understood. 

Of course, ultimately, the full signal is given by the 
contributions from all the relevant degrees of freedom. 


VI. OUTLOOK 

As the above examples have shown, the method pre¬ 
sented here is very powerful both in scope and predic¬ 
tive potential. We have, for example, explicitly shown 
that various hidden orders could be unambiguously iden¬ 
tified. Moreover, as we tried to emphasize, this approach 
offers the advantage of possibly helping with unbiased 
data analysis since all possible contributions to the RIXS 
signal can in principle be systematically enumerated. 

With this theory in hand, where should one look next? 
As proposed here, NiGa 2 S 4 of course appears as a natural 
material to investigate with RIXS or REXS. In particu¬ 
lar, thanks to S' = 1 on Ni 2+ one expects “direct RIXS” 
processes to be involved and therefore a strong signal. 
The current resolution on RIXS instruments —of about 
130 meV— is too low to detect a sizable signal-to-noise 
ratio for a material where the exchange has been esti¬ 
mated to lie at around J ~ 7 meV (as boldly estimated 
from a Curie-Weiss temperature of |@cwl ~ 80 K [15]). 
However, since static order is expected (at higher tern- 


peratures) [22] , Bragg peaks should appear in REXS (see 
Fig#))- Spin chain materials like LiCuV 04 and others 
m, while perhaps even more promising in terms of con¬ 
fidence in the realization of a nematic state, will have to 
await the next generation of RIXS instruments, as their 
exchange energies are also relatively low (~ 30 K). Per¬ 
haps, at this point, high-quality data (like in the cuprates 
and iridates) would be worth re-investigating in light of 
all the possibilities which our work unearthed. One can, 
for example, imagine looking for signs of some of the 
“stranger” correlation functions presented in Table[I] An¬ 
other exciting direction, briefly mentioned in Section [Vj 
is that of pnictide materials, as RIXS may help contribute 
to the effort of pinning down the origin of the observed 
nematic order. Finally, most electrifying would perhaps 
be the detection of chiral order in putative spin liquids on 
the kagome lattice [4TH43] or the possible appearance of 
spin quadrupolar correlations (in the absence of dipolar 
ones) in L^-aBa^C^O^ should it display features of a 
spin density wave glass @3]. 

With RIXS taking the central stage in various classes 
of systems, and new resolution-improved machines on the 
horizon, the future seems bright for refining our under¬ 


standing of and discovering yet new physics in complex 
materials amenable to RIXS. And with these general re¬ 
sults and derivation in this broad setting, we hope to 
guide experiments as well as theory in this endeavor. It 
is also our hope to have somewhat demystified the un¬ 
derstanding of RIXS for non-experts of microscopic cal¬ 
culations. 
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Appendix A: Electromagnetic field 

As mentioned in the main text, the electromagnetic 
vector potential at point r may be expanded in plane 
waves 

AW = E/^EK'-*' + t,.), (AD 

k e 
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where H is Planck’s constant, eo is the vacuum dielec¬ 
tric polarization, Wk = k = c|k|, with c the speed of 
light, V is the volume in which the electromagnetic field 
is confined, e has nonzero components only along (real) 
vectors perpendicular to k, and we define Ak. e (r) = 
£ *°ke e_?k r + h.c.. a ke is the creation operator of a 
photon of momentum k and polarization (helicity) e. 
This “expansion” introduces (and defines) the polariza¬ 
tion vector £ which encodes the vectorial (in the sense of 
a tensor of rank one) nature of the 5=1 field A. We 
return to the symmetry transformation rules of A and e 
in Appendix |Bj 


Appendix B: Transformation rules 


• zero spins: 

1 £ ' x 1 £ = (V £ + 1 £ ' X£ + ; (Cl) 

• one spin: 

1 £ ' xl E x I s - = (V £ + 1 £ ' X£ + 2^) x l Si ; (C2) 

• two spins: 

1 £ ' xl e x I s ’ x I s * = (V £ + 1 £ ' X£ + 2l £ '’ £ l) (C3) 

x (o S ‘' S i -J- fSixSj + 2 [s,,s j A 


The (pseudo-)vector of spin operators S r transforms 
under a spatial operation R and time reversal (TR) ac¬ 
cording to 


R : S r 
TR : S r 


det (R)R ■ S#., 
-S P 


(Bl) 


regardless of the value of 5(5 + 1) [45]. The vector po¬ 
tential A transforms as 


R : A(r) -x R ■ A (R ■ r) 

TR : A(r) -X -A(r) 


(B2) 


so that the polarization £ transforms according to 


three spins: 


1 £ ' x 1 £ x I s - x l s ^ x l Sfc 


(C4) 


(0 £ '- £ + l £ ' X£ + 2l £ '’ £ l) 

x (o SlSj + 1 s - xf h + 2l s - s jlj x l Sfc 

(0 £ '- £ + l £ ' X£ + 2l £ '’ £ l) 

x ^(Si-SjJSfc + 0 (S i xS J ')-S fc + 1 (S,xS,)xS t 

+ 2 [(S i xS J ),S fc ] + X [S* ,Sj ]-Sfe + 2 [Si ,Sj] x Sfc 
_j_ 3 [[Si ,Sj-] ,SfcI 5 


R : £ —X R • £ 

TR : £ -X —£* 


(B3) 


Note that the definition of the polarization sometimes 
differs by, e.g., a factor of i, and the polarization is then 
“even” (times complex conjugation) under the time re¬ 
versal operation. 

If the spatial symmetry group contains all spherical 
operations (which contain in particular all 50(3) oper¬ 
ations), S and £ transform under the “L = 1” repre¬ 
sentation of 50(3) (regardless of the value of 5(5 + 1)). 
Note that here, the name “L” is purely formal. Using 
the notation from Ref. nn for the full rotational symme¬ 
try group 11 D” (50(3) C D), e and S transform under 
the Oj” and D+ representations, respectively, where ± 
indicate parity under the inversion transformation. 


Appendix C: Derivation of Table I 


where the definition of the double brackets has been ex¬ 
tended to: 

! ([u, v] • w) p = £ Ju, v] F w„ 

([u, v] X w) Mp = A UAp[u, v]^U) A (C5) 

([[u,v],w])^ A = [u, v] pl/ UlA 

Only products of terms belonging to the same represen¬ 
tation will have a contribution in the “final” 0 represen¬ 
tation (by contracting all the indices). 

Explicitly, the operator obtained for all the terms in 
Table U reads: 

T = (s'* • e) KrSi • Sj + 00 , 2 ( 8 , x Sj) • S fc ] (C 6 ) 
+ (s'* x e) • [aipS, + ai^Sj x Sj + ai , 3 (Si • Sj) S fc 
+ai , 4 (Si x Sj) x S fc + ai )5 [Si, Sj] ■ S fe ] 

+ [£'*,£] ( 02,1 [Si, SJ + a 2 , 2 [Si x Sj, S fc ] 

+a 2 ,3 [Si,Sj] x S fc ). 


In the equations below, the numbers are representation 
labels (L = 0,1,2,... associated to D±), and the super¬ 
scripts schematically show basis elements (in the form of 
tensors) in terms of the original terms in the products. 
Products of representations for 


Appendix D: Lower symmetry 

It has been pointed out [14] that, even when spin-orbit 
coupling is negligible in the low-energy manifold, spin- 
orbit is always very strong in core levels, and may lead 
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to anisotropies in the RIXS signal. The derivation pro¬ 
vided in the main text is readily generalized to the case 
of discrete “spin” symmetries. 

With the help of the tables found in Refs. OEJ I35H351 
one may build bases for the representations, generaliz¬ 
ing those of rotationally-invariant systems. The formula 
which generalizes Eq. ( |C6[ ) is: 

T n = Y,Y, a rJ £T ' 1 - Sr ’ 1 ' ( D1 ) 

r i 

where the sum proceeds over all irreducible representa¬ 
tions T of the point symmetry group at site R (that of 
the core hole), l indexes the multiplicity of the represen¬ 
tation T, and the dot product represents a symmetric 
contraction of all indices. 


Appendix E: Higher multipoles 

As mentioned in the main text, the Hamiltonian at 
first order in A is actually 

■ a + $ f ^ T 4> r ■ (V x A r ) . 

m 2 m 

(El) 

In the main text, only the first term was considered. In 
the spirit of the derivation provided in the main text, 
where experimental parameters are kept explicit, to treat 
the second term, one should consider the couplings to 
k x e and k' x e', much as we did to e and e’ upon con¬ 
sidering the term linear in p. One should also note that 
“higher multipoles” will also arise from the expansion of 
the exponential, e lk Sr = 1 + «k • dr — |(k • (Sr) 2 + • • •. 


h ' = E 


Appendix F: Details of the cross-section derivations 
1. Spin nematic in S = 1 triangular magnets 


Furthermore, 

(Si • S,) 2 = —Si ■ Sj + \ 'E,{S?,Sr}{S?,S]}, (F4) 

fl,is 


and 

(F5) 

A 4 ) 1 '' 

= £(^W) 2 + \ 

with 


'{S*,Sy}=alp r +pta r 

= -*(&-$) 
{ S y,S*} = -i(-a r + at) 
(S?) 2 = 1 - Ptfc 
(Sj!) 2 = 1 — ala r 
,(4 Z ) 2 = c&a I +ftp T 


(F6) 


Using the rules Eq. ( |F3| ) and then keeping only terms 
quadratic in the boson operators ay, aj, j3 r and (i.e. 
neglecting interactions between the bosons), we arrive at 

H = \ (•Jl - J 2) Y EE l'lr T ll+R r , + VrVr+R„ 
ri—a,0 r n 

+y E EE[^+r„+m! + r „ 

?7=q,/ 3 r n 

-y E E E [nivr + ril +R Jlr+R„ 

r)—a,(3 r n 

= y ( J 1 - J i) Y Y [4^-k + VkTi-k 

r)=Oi,fl k 

+(Ji7k - 6J 2 ) EE vb)k, (F7) 

r)—a.,(3 k 


Following Ref. [19] (the calculation is performed there 
in the antiferroquadrupolar phase), we introduce the 
bosonic operators a r and /S r , and the Fock space vacuum 
such that 

f I'S'r = 0) = | vac) 

\ l-S'r = ±!> = 75(4 ± * 4 ) I vac) 

and 

f S* = aj + a r 

\ S V r = 4 + 4 , (F2) 

[s; = -i{alPr - fiar) 

with the constraint that there should be no more than 
one boson per site. This in particular implies, in real 
space: 

i a T = 4 = a r/3 r = 4 a r = 0 
a T pt = p r al = 0 . (F3) 

a r al = Ml = 1 - 4« r - 44 


where 7 k = 2(cos k x + cos {\{k x + y/3k y )) + cos {\{k x — 
\/3 k.y))). With the Bogoliubov transformation r = 
coshyPk + sinhyp_k, we obtain 

H = Y E Wk hp k, (F8) 

P=P°‘,P f> k 


where we have defined: 

Wk = 2 [A k cosh 2y + Rk sinh 2£ k ] = 2 A 2 . - R 2 , (F9) 


with 

f-Ak = 7k — 6 J 2 ) 

[Rk = f (4 - Ji) 

if 

J Ak sinh 2£ k + Rk cosh 2<y = 0 
1 (A k cosh 2y + Rk sinh 2£ k ) 2 = A\ - R£ 


(F10) 


(Fll) 
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This yields: 


sinh 2 2a = a^bI 
cosh 2 2a = A 2 A * B 2 


Since Wx coshx > 0, 


cosh 2a = 


K 


(F12) 




K-Bl' Sinh2Ck_ 

(F13) 


The transition operator Eq. 
Fourier space: 


takes the form, in 


• in the nematic phase: 

X {(s + s+)-„,~ k (s + s+)„, k ) 

0(w 2 — v+(k — £7r) 2 ) 


(F22) 


^E 

e=±l 


— v+(k — eir) 2 


2-1/K+ 


+^E e 


uj 2 — v+ ( k — 7r(e- — e'M) 


x \l oj 2 — v+ ( k — 7r(e^ — e'M) 


-K+ + 1/K+-2 


j( s - u ,,~k s ^,k'> = gapped (F23) 

I^-u-kXu.k) (x w _|_ v+k) + S(co — v+k)) (F24) 


T k = K ( 0) S(k) - EE [ pI(A m cosh a + Bf, sinh a) 

M p=p(t 

+/3-k(-^A‘ sinh a + Bfj, cosh a) , (F14) 


I^-v.-kSZ.k) (x w (<5( w + u + /c) + S(lu — v + k)) (F25) 

y- 0(w 2 - v+(k - e7r(| - M)) 2 ) 

+ 2^ /-2-AT+ 

e=±1 Jw 2 - v 2 (/c — £7r(i — M)) 2 


where 

(F15) 
(F16) 

(F17) 

(F18) 
(F19) 
(F20) 
(F21) 

Plugging this into the expression for the cross section: 

£ 2(J 

oc y 1cosh £ k + B M sinh £ k | £(w-u; k ), 

y=x,z 

we arrive at the result given in the main text. 


« (0) = a 0 (e ■ e'*) 


cl 1 ) = 


^1 ^-fl(7T^ & 


e 1 - 


012 

- 

3 <W 

A 2 ) 

^xy 

+ 

in z 

A 2 ) 

n yz 

~ 

in x 

A 2 ) 

^xy 

- 

in z 

A 2 ) 

^yz 

+ 

IK X - 


2. Vector chirality and bond nematic in S = 1/2 
Ji — J 2 chains 


Equal-time and real-space correlation functions are 
given in Refs. [351 1313 ■ Here, we find the following con¬ 
tributions to the cross section: 


• in the vector chiral phase: 

X(x^,- k xZ, k ) (F26) 

oc Atv 3 ( 5(ui + v+k) + S(oj — v+k)) 

+B J2 [©(w 2 -u 2 (fc-e27rM-e , 7r) 2 )] 

e,e'=±l 

/-4A-+-2 

xww 2 - v+(k — e2irM — e'le) 2 

X( s -*,-k s *,k) (x w (5{u> + v+k) + S(oj — v+k)) (F27) 

0(w 2 — t; 2 (k — eQ) 2 ) 

-2-l/(4A~a 

£=±1 Juj 2 — v+(k — eQ) 2 


X< s -«,-k s Z,k) (x ^ 


j{(s+s + )_„ ; _ fc (s , +s ,+ )„ |Je ) (F28) 

0(w 2 — v+(k — e2Q) 2 ) 


E 

e=± 


1 \Ju 2 — v+(k — e2Q) 


- 2 - 1 /K+ ' 


In all the above, A and B are constants, and Q = § — 

<W- 

Note, in particular, that, since K increases mono- 
tonically between 1/2 and 1 ( K(M = 0) = 1/2 and 
K(M = 1/2) = 1), and K+ = K( 1 + K%), K+ > 1/2. 
Moreover, the bosonization approach is valid only “not 
too close” from the saturation limit M = 1/2, and in the 
weak coupling regime v ~ J 2 . So, in particular: 


2 - Y > 3/2 
K+ + + Y - 2 > 0 
2 - I<+ > 0 
4 K+ - 2 > 0 


• r T+ ^ -l + y/l+8Jl/(irv) 

if A < 7 TV - 271- ■ 


2 - 


1 > 15 

4A+ — 8 


(F29) 

Note that K < nv 1+ -\/ 1 ^ Jl /(^-») ^ a j wa y g ^ rue f or A < 

1. 
































